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BOtKDAKY  STABILIZATKH  OF  THIN  ELASTIC  PATES 

John  E.  Lagnese 

_  AFOSRTR.  8  7-  1  562 

Department  of  Mathematics  J  ** 

Georgetown  University 
*  Washington,  DC  20057 


1.  Introduction.  In  this  paper  we  shall  consider  the  question  of 
uniform  stabilization  of  thin,  elastic  plates  through  the  action  of  forces 
and  moments  on  the  edge  of  the  plate  (  or  on  a  part  of  the  edge  of  the 
plate).  Two  particular  plate  models  will  be  considered:  The  classical 
fourth  order  Kt rchhoff  model ,  but  incorporating  rotational  inertia,  and  the 
"sixth  order"  Hindi in*Timoshenho  model.  The  difference  in  the  two  models, 
from  a  physical  point  of  view,  is  that  the  M-T  model  incorporates 
transuerse  shear  effects  while  the  Kirchhoff  model  does  not.  Actually,  the 
M-T  model  is  a  hyperbolic  system  of  three  coupled  second  order  partial 

differential  equations  in  two  independent  variables.  The  unknowns,  denoted 

<<V  v 

by  w.  <?  are  the  vertical  component  w  of  displacement  and  angles  '//.  <f 

which  are  measures  of  the  amount  of  transverse  shear.  The  three  equations 
are  coupled  through  terms  which  are  multiples  of  a  factor  K  called  the 
coefficient  of  elasticity  In  shear^It  is  well  known  that,  as  K-**>,  the 
limit  of  w=w^  is  formally  a  solution\of  the  fourth  order  Kirchhoff  plate 
equation  (2.1)  below  (see  Remark  2.1  telow.  The  manner  in  which  certain 

conservative  M-T  systems  converge  to  tlip  corresponding  Kirchhoff  systems 

\ 

has  been  made  precise  in  J.  LAGNESE-J.L.\  LIONS  [10,  Chapter  II].) 

For  the  M-T  model,  it  will  be  showri  that  one  can  find  simple,  natural 
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feedback  laws  for  the  edge  forces  and  moments  which  lead  to  an  explicit, 
uniform  decay  rate  u=u(K)  for  solutions  to  the  closed-loop  system  (under 
some  geometric  conditions  on  the  shape  of  the  plate).  It  will  be  seen, 
however,  that  u( K)-«0  as  K-*».  On  the  other  hand,  we  can  show  that,  as  Kh», 
the  limit  of  w=w^  exists  and  Is  the  solution  to  a  dlsstpatlue  closed-loop 
Kirchhoff  system.  One  may  then  prove  directly  that  the  limit  system 

|f 

likewise  has  a  uniform  decay  rate  u  .  Thus  we  have  the  puzzling  situation 

|f 

that  the  limit  of  the  decay  rate  u(K)  is  not  the  decay  rate  u  for  the 
limit  system. 

A  related  paradox  was  encountered  in  [10]  in  connection  with  exact 
boundary  controllability  of  the  M-T  model.  There  it  was  shown  that 
solutions  of  the  M-T  model  are  exactly  controllable  in  a  finite  time 
Tq=Tq(K)  .  where  Tq(K)hE0  as  K-<“.  Nevertheless,  the  limit  of  the  M-T  system 
is  a  Kirchhoff  system  which  can  be  shown  to  be  exactly  controllable  in  a 

|f 

finite  time  T^.  It  should  be  remarked,  however,  that  these  anomolies  may  be 
a  result  of  the  methods  employed  rather  that  anything  intrinsic  to  the 
models  themselves. 

For  related  work  on  uniform  boundary  stabilization  of  (mainly) 
hyperbolic  equations  we  refer  to  G.  CHEN  [1.2],  V.  K0M0RNIK-E.  ZUAZUA  [6], 
J.  LAGNESE  [7,8.9],  I.  LASIECXA-R  TRIGGIANI  [11],  and  J.L.  LIONS  [12].  For 
results  on  stabilization  of  beam  problems  with  end  damping,  see  G.  CHEN  et 
al  [3.4]  and  J.U.  KIM  and  Y.  RENARDY  [5]. 

2.  Models  to  be  Considered.  Consider  a  homogeneous,  isotropic,  thin 
plate  of  uniform  thickness  h.  Points  within  the  body  will  be  represented  by 
rectangular  coordinates  (x.y.w).  It  is  assumed  that  the  plate  has  a  middle 
surface  midway  between  its  faces  which,  in  equilibrium,  occupies  the  region 


K 


2 


Q  of  the  plane  w=0.  We  denote  by  w(x.y.t)  the  v-component  of  the 
displacement  vector  at  time  t  of  the  point  which,  when  the  plate  Is  In 
equilibrium,  has  coordinates  (x,y,0). 

It  is  customary  in  thin  plate  theory  to  assume  that  the  transverse 
normal  stress  a is  negllgable  compared  to  the  other  stresses.  It  is 
further  assumed  that  all  components  of  displacement  are  "sufficiently 
small"  to  Justify  linearization.  (The  reader  may  consult  any  number  of 
works  for  details  of  the  linear  theory  of  elastic  plates.  A  particularly 
accessible  source  is  the  book  of  K  WASHIZU  [13].  Heuristic  derivations  of  a 
number  of  plate  models,  both  linear  and  nonlinear,  may  be  found  in  Chapter 
I  of  [10].) 

The  Kl rchhoff  model  results  if  one  further  assumes  that  the  linear 
filiments  of  the  plate  initially  perpendicualr  to  the  middle  surface  remain 
straight  and  perpendicular  to  the  deformed  middle  surface  and  undergo 
neither  contraction  nor  extension.  (Thus  transuerse  shear  is  neglected.) 
These  assumptions  lead  to  the  following  partial  differential  equation  for 
w(x.y. t): 

(2.1)  phw"  -  Aw"  +  DA2w  =0  in  0x(0.»). 

d 

where  '  =  g-p  A  is  the  ordinary  Laplacian  in  (x.y)  variables,  p  is  the 

3  2 

surface  density  per  unit  area  and  D=Eh  /12(l-p  )  is  called  the  modulus  of 
flexural  rigidity  (here  assumed  to  be  constant,  p  is  Poisson’s  ratio  and 
satisfies  0<p<l/2,  and  E  is  Young’s  modulus. 

We  assume  that  the  (piecewi se-smooth)  boundary  T=d0  consists  of  two 
disjoint  parts:  r=roUTj  with  and  Tj  relatively  open  in  T.  We  denote  by 

i>  the  unit  normal  to  T  pointing  towards  the  exterior  of  Q,  and  we  set 
^n=1~nx  (O,00),  S^r.x^,10) .  The  boundary  conditions  to  be  considered  are 


(2.2)  «  =  §“  -  0  on  20 

(2.3)  D[Aw  ♦  (l-ji)Bjir]  =  Vj  on  2j 

(2.4)  d[1±i  ♦  o-m)b,»]  -  mr  ‘ 


on  20. 


12  dv 


where 


no  2  2 

B.w  =  2uiu0w  -  u.w  -  u0w 

1  1  2  xy  lyy  2  xx 

B0w  =  -  i>o)w  ♦  ~  w  )] 

2  5tlv  1  2'  xy  1  2V  yy  xx' J 

and  where  Tsf-Dg.Uj)  is  a  unit  tangent  vector.  We  have  used  (and  will 

continue  to  use)  the  notation  w  =3w/9x,  w  =o  w/dxdy,  etc. 

x  xy 

The  boundary  conditions  (2.2)  are  the  so-called  geometric  boundary 
conditions  and  mean  that  the  plate  is  clamped  along  Tq.  Conditions  (2.3) 
and  (2.4)  are  the  natural,  or  mechanical  boundary  conditions .  Condition 
(2.3)  corresponds  to  bending  about  the  axis  formed  by  the  tangent  to  T. 
while  (2.4)  incorporates  both  shear  force  in  the  w  direction  and  twisting 
around  the  normal  to  T.  v^  and  are  the  control  uariables  through  which 
the  system  is  to  be  stabilized.  Of  course,  to  uniquely  determine  the 
motion,  initial  conditions 
(2.5)  w(0)  =  w°.  w' (0)  =  w1  in  0 

must  also  be  specified. 

The  Mindlin-Timoshenko  model  arises  if  the  Kirchhoff  hypothesis  is 
weakened  by  removing  the  assumption  that  the  filiments  of  the  plate  remain 
perpendicular  to  the  deformed  middle  surface,  but  retain  the  assumption 
that  these  filiments  remain  straight  and  undergo  no  strain  in  deformation. 
It  is  a  system  for  three  functions  w.  >p,  < p  where  w  has  the  same  meaning  as 

fV 

before,  ^  ~  -(w  +4'),  <p  ~  - (w  +tf>)  and  where  4>.  <p  represent  angles  of 
x  y 

rotation  of  the  cross  sections  x=const.,  y=const.  containing  the  filiment 
which,  when  the  plate  is  In  equilibrium,  is  perpendicular  to  the  middle 
surface  at  the  point  (x.y.O).  The  system  is  the  following: 


4 


(2.6) 


i 

1 


a 


(2.6)  +"  -  ^  ^  *  K(*  *  wx>  *  °- 

(2.7)  £|-  *"  -  D^*  ^  ^  +xyl  4  K<»  4  ",>  -  0. 

(2.8)  phw"  -  K[(^wx)x  ♦  (^wy)y]  =  0 

in  0x(0.“).  The  geome  trie  boundary  conditions  for  (2.6)-(2.8)  are 


(2.9) 


^/  =  <?  =  w  =  0  on  Iq, 


while  the  mechanical  boundary  conditions .  applied  on  2j .  are 

(2.10)  Dt>i*x  +  +  ^x^2^  =  VT 

(2.11)  9[  v2lf>y  +  **v2*x  +  ^x^l^  =  v2‘ 

(2.12)  K[Wi)  +  +  i>2*]  =  v3 

where  usfOj.Dg).  Conditions  (2.10)  and  (2.11)  correspond  to  twisting 
moments  while  (2.12)  corresponds  to  shear  force  in  the  w-direction.  The 
Vj’s  are  the  control  variables. 

The  Initial  data  for  (2.6)-(2.12)  are 

*(0)  =  *°.  *‘(0)  = 

(2.13)  <p(0)  *  <f°,  <?'  (0)  =  (p1 , 

w(0)  =  w°.  w ' (0)  =  w1  in  0. 


The  constant  K  which  appears  in  the  above  equations  satisfies  K/h  = 
kE/2(l+p).  K/h  is  called  the  coefficient  of  elasticity  In  shear;  k>0  is  the 
shear  correction  coefficient .  In  Section  4  we  shall  study  the  convergence 
of  the  system  (2 . 6)  — (2 . 13)  as  K-*0. 

Remark  2.1.  Equations  (2.6)-(2.8)  may  be  formally  uncoupled  (see. 
e.g. .  [10.  Chapter  I])  to  yield  the  following  equation  for  the  vertical 
displacement  w: 

3  3 

(2.14)  phw"  -  Aw"  +  DA2w  +  jp|j2  w""  -  DAw"J  =  0. 

Thus  (2.1)  is  the  formal  limit  of  (2.14)  as  K-*». 

3.  Uniform  Stabilization  of  the  Hindi  in-Timoshenko  Model.  Let  4>.  <p , 
w  be  a  solution  to  (2.6)-(2. 12) .  The  strain  energy  at  time  t  corresponding 


5 


to  this  solution  Is  given  by 

(3.1)  y(t)  =  \  jr  to2  +  i?  ♦  2MVy  +  *x)2J 

♦  K[(*  +  wx)2  +  (<f  *  wy)2]}dxdy, 

and  the  kinetic  energy  is 

2 

(3.2)  *(t)=  &  s  [W2  +  Y5<^'2+  f  2)]dxdy . 

z  n  id. 

The  total  energy  E(t)  is  defined  to  be 
E(t)  =  <f{  t)  +  *(t). 

We  wish  to  determine  feedback  controls  Vj  in  (2. 10)-(2. 12)  such  that  the 
resulting  closed-loop  system  is  uniformly  asymptotically  stable,  i.e., 

(3.3)  E(t)  *  Ce^^CO) 
for  some  posltlue  constant  u. 

In  order  to  determine  candidates  for  stabilizing  feedbacks  we 

differentiate  E(t)  in  t.  We  obtain 

E'(t)  =  S  (Vj*'  +  v2<p‘  +  v3w')dT. 
fl 

Therefore  the  feedbacks 

(3.4)  Vj  =  -kj+’ .  »2  =  -k2f  '.  v3  =  -k3»’ 

with  k^O  introduce  dissipation  into  the  system  in  the  sense  that  E(t)  is 

00 

nonincreasing.  Actually,  we  can  prove  that  if  k^L  (Tj)  satisfies 
kj(*.y)2kQ>0  a.e.  on  Tj .  1=1. 2. 3.  then  under  some  geometric  conditions  on  T 
we  do.  in  fact,  have  (3.3).  However,  the  required  geometric  conditions 
force  if  T  is  smooth  and  thus  exclude  simply  connected  regions  Q 

with  smooth  boundaries  if  1=1,2,  (which  we  explicitly  assumed. 

However,  see  Remark  3.3  below  concerning  the  assumed  regularity  of 
solutions . ) 

In  order  to  weaken  the  geometric  restrictions  and,  simultaneously,  to 
obtain  explicit  expressions  for  C,  u  (in  terms  of  the  parameters  of  the 


system  (2.6)-(2.8)  and  the  geometry  of  T).  we  alter  the  feedbacks  (3.4)  in 
a  way  which  Incorporates  the  geometric  features  of  T.  (This  Idea  was  first 
Introduced  in  [6]  in  connection  with  stabilization  of  solutions  of  wave 
equations.)  To  this  end,  let  us  set 
m(X)  =  X  -  XQ 

2 

where  X=(x,y),  Xo=^xO’yO^’  XQ  a  fixed  point  of  R  .  We  assume  that  XQ  may  be 
chosen  such  that 

(3.5)  m*u  £  0  on  Tj . 

(3.6)  m*u  £  0  on  Tq. 

We  then  define  feedbacks  by  setting 

(3.7)  v  j  =  -(m*u)kj*',  v2  =  -(nfu)^1 .  v3  =  -(m*u)k3w1 
where  k  CL^Tj).  ki(x,y)^kQ>0  a.e.  on  .  1=1, 2. 3.  We  then  have 
THEOREM  1.  Assume  that  Tq,  Tj  satisfy  (3.5),  (3.6)  and  that  the  controls 
Vj  tn  (2. 10)-(2. 11)  are  given  by  (3.7).  Let  the  initial  data  (2.13)  satisfy 

+°,  <f°.  w°c  H!(n).  *°=  v>°=  w°=  o  on  r0. 

+1.  f1.  w*€  L2(0). 

Then  there  exists  a  constant  u=u(K)  (luhose  explicit form  is  giuen  in  (3.8), 

(3.11)  below)  such  that 
00 

S  E(s)ds  s  r-  E(0). 

0  w 
00  00 

S  E(s)ds  1  e^V  E(s)ds ,  t  ^  0. 

t  0 

COROLLARY.  Under  the  conditions  of  Theorem  1,  me  haue 
E(t)  *  e"Wt+1,  t  *  1/u. 

Remark.  3.1.  The  controls  v^  of  (3.7)  satisfy  v^€  L2(2j),  and  the 

problem  (2.6)-(2. 13) ,  (3.7)  is  well-posed  in  the  space  (hI  (n))^x(L2(D) )^, 

0 


where 


Proof  of  Theorem  1.  The  proof  Is  very  lengthy  and  technical  and  will 
only  be  sketched  here.  Complete  details  will  be  published  elsewhere. 

For  0<e<l,  define 

Fe(t)  =  Yj(t)  +  (l-e)Y2(t)  4  tY3(t) 


where 


Yj(t)  =  ph 


■  2  2 

(t)  ,m*W(t))  +  (t).m*v^(t))  +  (w' (t),m*vw(t)) 


y2(0  »  f|-[(*,(t).*(t))  +  (-f'(t).^(t)) 
Y3(t)  =  ph(w  (t).w(t)) 


and  where 


=  /  v>(x.y)^(x.y)dxdy. 

n 


Also,  for  6>0.  define 


E£  fi(t)  =  E(t)  +  6Ffc(t). 

Since  | ( t ) |£CE(t)  with  C  independent  of  a  and  K  (provided  1QK0>0),  we 
have  Efc  ^(t)^0  provided  6C<1.  One  can  prove  the  following  (and  this  Is  the 
key  point)'-  There  are  positive  constants  c,  Xq  and  X^  (which  will  be  made 
explicit  below)  such  that,  if  we  choose 


(3-8) 


0  <  e  < 


1 

2(14KX0)- 


0  <  6  < 


2c  c 

eph4cXj ’ 


then 

(3.9)  grEfc  fi(t)  *  -6eE(t).  t*0. 

The  constants  c.  Xq,  X^  are  defined  as  follows. 

c  =  max  sup  k. (x.y) . 

1  F1 

In  order  to  define  Xq  and  Xj .  we  introduce  the  notation  (which  will  also  be 
used  in  Section  4) 

=  Dj  +  +  ,^(^y+  fx)2  +  2M+X<f>y]dxdy. 


•*-  -  ^  X.  ‘A-  -  •A-  J  A  “* 


.Va  ’ji 


'_*»  *  *  *  •  * 


8 


«Y  (*.*.*)  «  S  ■*1»{D[^X  ♦  * y  ♦  ~2^*y+  *X)  +  + 

1  rl 

k[(*  +  wx)2  +  (*  ♦  wy)2]}dr. 

b(^.<#>,w)  =  S  m'udt^2  +  k^-p2  +  k3»2)dr. 

F1 

Then  Xq  and  X^  are  defined  by  the  inequalities  (which  incorporate  the 
geometry  of  0) 

H'P2  +  <f2)dxdy  i  X_an(^.<f).  V^.^€h!  ((7) 

Q  U  U  *0 

(this  last  inequality  is  not  true  if  Tq=$) . 

bfnrv^.m^v^.m'Vw)  i  X.  [a,.^.^)  +  a_,  (^.f.w)],  ,  <p ,  w€hI  (f7)nH^  (f7)  - 

l  u  i  j  iQ 

The  conclusions  of  Theorem  1  are  obtained  from  (3.9)  as  follows.  Let 

(5>0.  multiply  (3.9)  by  e  and  integrate  in  t  from  t  to  °°.  After  an 

integration  by  parts  we  obtain  (since  E  r(t)£0) 

£  » o 


(3.10) 


Define 


(3.11) 


r  -Psw  ,  \  ~Ptr  ,  _-P  t  1+6C 
|eME(s)ds^e  Efe.6*e  ~ 


E(t). 


and  let  piO  in  (3.10).  The  result  is 


S  E(s)ds  *  ^E(t).  t  *  0. 

The  last  estimate  implies  the  conclusions  of  Theorem  1. 

The  Corollary  to  Theorem  1  is  a  consequence  of  the  inequality 
00 

(3.12)  tE(  t+T)  i  S  E(s)ds  £  -e^^fO).  Vr>0. 

t  W 

In  fact,  from  (3.12)  we  obtain 

£{  t+T )  s  W”<t4T>. 

'  ’  urr 

The  factor  e^/urr  has  its  minimum  at  t=1/w  and,  for  this  value  of  r,  we 
obtain  the  estimate  of  the  Corollary. 

Remark  3.2.  From  (3.8)  and  (3.11)  it  follows  that  w=c>(K)-0  as  K-*1 


Remark  3.3.  The  proof  of  Theorem  1  requires  that  solutions  be 
classical .  1 .e.  H^(0) .  However,  there  Is  a  difficulty  due  to  the 
possibility  of  singularities  at  points  of  j.  which  may  preclude  the 

existence  of  classical  solutions.  Such  singularities  can  occur  regardless 
of  the  regularity  of  the  initial  data.  Thus,  while  Theorem  1  is  certainly 
valid  if  it  is  proved  in  the  opposite  case  only  if  one  has  a 

priori  knowledge  of  the  assumed  regularity. 

4L  Limiting  Behavior  as  K  -*■.  We  introduce  H=L^(fi)  with  the 
standard  scalar  product,  and 

V  =  .pCH^fi).  *  =  0  on  rQ} . 

We  further  define 


a  0(*. +  (*>„.<#>„)  +  + 


x  x 


y  y  x  y' 

1-i 


+  VV  *x)]' 


AAA 


a^.f  .w;^.<f  .w)  =  (^wv.Nk+wv)  +  .  <P+wv  ) , 


AAA 


=  J  +  k^W  +  k^ww^F. 

ri 

2  2 

+  +  (»■*) 


c('4/.<f>.*;'P.>p.w)  =  ph 


The  solution  to  the  boundary  value  problem  (2.6)-(2. 13) ,  (3.7)  may  be 
defined  via  the  variational  formulation 


b(^’  ,if>‘  ,w’  w)  =  0.  V^.f.  w€V. 


(4.1) 

d  A  A 

dT**’ 

(4.2) 

*.*,w€C([0.T];V). 

(4  3) 

[(m*n)k  ]^'  |_  , 

1 

(4.4) 

*(0)  =  *°.  *(0)  = 

♦  '(0)  =  f1.  sp- (0) 

If  we  assume  that 


It  way  be  shown  that  the  problem  (f  l)-(f  <)  is  well  set  We  better  i  •  % 


solution  by  4^.  w^  in  order  to  emphasize  the  dependence  or.  K  Wr  nor 
wish  to  study  the  behavior  of  w^  as  K-«*  To  do  so  wr  adap'  *r.» 

technique  of  [10.  Chapter  II]  to  the  present  problem 

AAA 

We  choose  <f>.  w  equal  to  ^(t),  <#>^(t).  w^.(  t )  .  resprr  tiven 
(4.1)  and  integrate  in  t  from  0  to  t.  (This  Is  formal  but  mnv  for  exam; 
be  justified  by  using  a  Galerkin  approximation  to  ^ ^  )  The  res 

is 


(4.6)  c(^(t).^(t).w^(t))  +  aQ(^K(  t )  ,^K(  t ) )  4 

Kal(V'K(t).^K(t).w1<(t))  +  S  S  m-Dfkj^2  4  k.^K2  4  k^K2)d’n- 

=  c(^1.«f>1.w1)  +  a0(^°.^°)  4  Kaj  (^°.<f0.w°) 

In  order  to  get  something  out  of  (4  6)  we  ncike  the  assum; >t  i  on  (;• 
addition  to  (4.5))  that 
(4  7)  a^AA*0)  =  0. 

i.e.  . 

A  A  Hi  (0).  w°«  H?  (0). 

0  0 

,o  o  „  o  o  _ 

J,  4  W  =  0.  ^  4  i  =  0 , 

x  y 

where 

(0)  =  (w:  wtH2^).  w  =  wt  =  0  on  T Q } 

In  addition,  we  assume  that 
wJ€V. 

Remarh  U  1.  A  physical  motivation  for  assumption  (47)  is  as 
follows.  The  strain  energy  associated  with  a  solution  w  of  the  Kirchhoff 
model  (2. l)-(2.5)  is  given  by 

(4  8)  t)  =  \  JT  D[w^  4  w^  4  2Mwxxwyy  4  2(1-M)w^]dxdy 

It  is  seen  that  (4.8)  is  formally  obtained  by  setting  'f'=-wx  and  -w  in 


UWAWWWUIIW 


WWW 


the  expression  (3.1)  for  the  strain  energy  associated  with  the  N-T  nodel. 
Thus,  hypothesis  (4.7)  way  be  Interpreted  as  the  requirement  that  the 
tnttlal  strain  energy  be  consistent  with  the  Kl rchhoff  model. 

Assuming  (4.7).  we  obtain  from  (4.6)  that 

(4.9)  are  bounded  in  L^O.TiV), 

00 

^K'  ^K'  WK  are  ^oun<^e(^  in  ^  (O.T;H), 

[(m-ujkj]14^.  [(m-ujkg]^.  [(m-u)^]^ 

are  bounded  in  1/^(0, )) , 

(4.10)  vR(*k-  wKx^'  *Ky^ 

are  bounded  in  L*°(0,T;H). 

From  (4.9),  (4.10)  follows  that  we  also  have 

WKx'  WKy  are  k°un£*e<*  ln  L*(0,T;H). 

Therefore,  we  may  choose  a  subsequence,  still  denoted  by  </»  w  such 

K  K  K 

that 

(4.11)  ^K>  *>K.  wK  -»  <p.  w  weak**  in  L*(0.T;V). 

(4.12)  4^.  <f^.  Wj ’  *>'  •  w’  weak**  in  l"°(0,T;H). 

(4.13)  (^K+  w^),  (<#>K+  »Ky)  -»  0  strongly  in  L*(0.T;H). 

(4.14)  [(m'ujkj]**^  -*  [(m»i>)kj]1V  weakly  in  L^O.TiL2^) )  . 

(4.15)  [(m*u)k2]\^  -♦  [(m*i»)k2],V  weakly  in  L2(0,T ;L2(rj ) ) . 

(4.16)  [ (m* i> -*  [(m'ujk^^w'  weakly  in  L2(0,T;L^(rj )) . 

From  (4. 11) . (4. 13)  it  follows  that 

(4.17)  *  +  w  =  0.  f  +  w  =  0, 

x  y 

(4.18)  r  +  w;  =  o.  *  +  w*  =  o. 

x  y 

Therefore 

(4.19)  w€L°(0.T;H?  (0)).  w,GL°‘,(0.T;h1  (0)). 

0  J0 

We  also  obtain  from  the  limiting  process  that 
w(0)  =  w°,  w' (0)=  w1 . 


K 


Pi 


12 
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Next,  in  (4.1)  we  choose  wGC  (17)  and  4>.  w  according  to 

0 

A  A  A  A 

>h  ~  -  W  ,  =  ~W  . 

T  x  y 

We  obtain  the  equation 


(4.20) 


dt  C^K^K*WK:'Wx'-WylW)  +  ^V'W'V  + 


b^K‘<fK,wK:-Wx,“Wy,w^  =  °*  Vw€Hp  (n)- 


Using  (4. 11)— (4. 16)  we  may  pass  to  the  limit  (in  the  sense  of 
distributions)  in  (4.20).  Let  us  write  out  the  limit  equation,  using 
(4.17).  (4.18).  We  obtain 

3 

(4.21)  ^tph(w‘ .w)  +  ^-{vw1  ,vw)]  + 

A  A  A  A  A 

DJ  [w  w  +  w  w  +  p(w  w  +  w  w  )  +  2(l-u)w  w  Idxdy  + 

P  L  xx  xx  yyyy  xx  yy  yy  xx'  v  xy  xyJ 

S  m*y(kjWWrx  +  k^w  +  k3w  w)dr  =  °>  VwCH^  (fl). 


(4  22) 


r(0)  =  w°€H?  (0).  w' (0)  =  w*Ch1  (0). 

0  0 


But  (4.21).  (4.22).  (4.19)  has  a  unique  solution  and  are  exactly  the 
variational  formulation  of  the  boundary  value  problem 


(4.23) 

(4.24) 


»"  -  Aw"  ♦  DA2w  =  0  in  iTx(O.T). 


dw  n 
dv 


on  1Q. 


(4.25)  D[Aw  ♦  (I-mJBjW]  =  -(■•vHkjVjW^  ♦  kgVy*  on  5j . 

3 

(4.26)  D[§£L  4  (1-M)B2w]  -  |j"  =  (n»-u)k3w  - 

g^m-ut-kju^  ♦  k^jWy)]  on  ly 


(4.27) 


w(0)  =  w°.  W(0)  =  w1 


in  (7. 


Remark  4.2.  It  may  be  shown  by  methods  similar  to  those  of  Section  3 
that  the  problem  (4 . 23)-(4 . 27 )  is  uniformly  asyrnpt ot Icall y  stable.  Thus 

ft 

this  problem  has  a  positive  energy  decay  rate  u  which  is  not  the  limit  of 
w(K)  as  K-*°  (defined  in  (38),  (3.11)). 

Remark  4.3.  If  kj=k2=k*.  the  right  sides  of  (4  25).  (4.26)  reduce  to 


(4.28)  -(m*t>)k* 
and 

(4.29)  (m'lOltjW'  -  |^(m*t>)k*  ^  ] 

3 

respectively.  When  the  inertia  term  -(ph  /12)Aw"  in  (4.23)  (and  the 
corresponding  term  in  (4.26))  are  absent,  uniform  asymptotic  stability  was 
proved  in  [9]  using  the  feedbacks  defined  in  (4.28),  (4.29)  with  k*=0. 
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